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In this note we consider the homogeneous and isotropic cosmology in the finite-range gravity 
theory recently proposed by Babak and Grishchuk. In this scenario the universe undergoes late 
time accelerated expansion if both the massive gravitons present in the model are tachyons. We 
carry out the phase space analysis of the system and show that the late-time acceleration is an 
attractor of the model. 
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I. INTRODUCTION 



Recently there has been renewed interest in the old 
idea that gravitons could have a small mass §. The fresh 
motivation for massive gravity may be attributed to the 
high precision experiments on the detection of gravita- 
tional waves like LIGO and LISA which, in a near future, 
may result in answering the question about the mass of 
gravitons ||, ||, §■ Secondly it seems to be natural to look 
for an alternative to GR to incorporate the recently ob- 
served " late time accelerated expansion" of the universe. 
Third, the higher dimensional gravity theories seem to 
naturally mimic the properties of four dimensional mas- 
sive gravity and should be given serious thought |fl4) . 
Last but not the least it is more than desirable to incor- 
porate the possibility of inflation without opting for an 
ad hoc mechanism to realize it. The massive gravity may 
or may not address all these issues. 

Till very recently, the masses of gravitons were thought 
to come from Fierz-Pauli mass term added to Einstein 
Hilbert action. The model faces the well known mass dis- 
continuity problem called Van Dam-Vcltman-Zakharov 
discontinuity!, g @, 1 1 0, 0, [jj ||, @, ^ [l| 0.1t 
was very recently noticed by Babak and Grishchuk that 
this problem is very specific to the choice of Fierz-Pauli 
mass term[^9). In field theoretic formulation of Babak 
and Grishchuk [fl0|, f20f , the mass term has the form, 
hh^h 



k)2h , where is the nonlinear gravita- 
tional field propagating in the Minkowski space and "h" 
denotes the tress of W v . Linear massive gravity is dis- 
cussed in Rcfs.^l], ^2[ For a specific choice of pa- 
rameters fci and k)2'. k\ + fc 2 — 0, the local predictions of 
the model have finite difference with those of GR though 
the mass term in the Lagrangian smoothly vanishes for 
ki — > 0. Babak and Grishchuk have demonstrated that 
everywhere on the parameter plane fo^) excluding the 
Fierz- Pauli straight line( fc 2 = — k\ ), the massive grav- 



ity theory proposed by them is free from discontinuity 
problem The parameters fci and fc 2 are related to the 
masses of gravitons: a - mass of spin-2 graviton and (3 - 
mass of the spin-0 graviton. The predictions of Babak- 
Grishchuk model are dramatic: their model has smooth 
limit to GR for a, (3 — > 0, the inclusion of a small mass 
term removes the black hole event horizon and makes the 
cosmological evolution oscillatory. Further, it is remark- 
able that if both the gravitons in the model are sacri- 
ficed to tachyons, the finite range gravity exhibits late 
time accelerated expansion of the universe. The mira- 
cle is done by the presence of massive scalar graviton. 
In case, the mass of scalar graviton is zero, Babak and 
Grishchuk have noticed that cosmology based upon the 
finite range gravity is identical to GR cosmology inde- 
pendently of the mass of spin-2 gravitonJl9|. Hence the 
additional spin-0 graviton plays a central roll in finite- 
range gravity. For earlier references on field theoretic 
formulation of gravity see Re fsfpij, ^5|, ^6) and the work 
of Logunov and collaborators fl27| ^The non- linear bigrav- 
ity theory recently discussed by Damour et al. seems to 
share many interesting features of the finite range gravity 



mentioned above 28 
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Once the Fierz-Pauli choice for mass parameters is ex- 
cluded, the success of the theory relies on existence of 
an additional scalar graviton. This is reminiscent of the 
manner in which a massive scalar boson was introduced 
by Veltman in the electroweak theory fFfl| , and sugges- 
tion contained in Ref. Jl8| that gauge bosons and massive 
graviton must necessarily be interpreted as multi-spin ob- 
jects. In particular, it is anticipated that a finite-range 
gravity should endow graviton with spin 0, 1 and 2 com- 
ponents. However, in this work we shall strictly confine 
to the framework of Babak and Grishchuk. 



In this note we investigate the phase space behavior of 
cosmological evolution in finite range gravity with both 
the gravitons behaving like tachyons. We demonstrate 
that the accelerated expansion is a late time attractor of 
the model. 
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EVOLUTION EQUATIONS IN FINITE 
RANGE GRAVITY 



The gravitational field in the massive gravity is de- 
scribed by a non-linear tensor field propagating in 
the Minkowski space with metric 7 A " / and with a massive 
term added to the GR part of the field Lagrangian which 
is quadratic in h?" (covariant derivative is taken with re- 
spect to the metric 7 MI/ ). The effective field equations 
have the form Jill 



(1) 



where T^ v is the matter energy momentum tensor and 
the mass tensor M M „ is given by 



M, 



(2) 



where 



In a homogeneous 
and isotropic situation h^ v depends upon time alone and 
the effective metric acquires the form 



ds 2 



b 2 (t)dt 2 - a 2 (t) (dx 2 + dy 2 + dz 2 ) 



(3) 



where a(t) and b(t) are given by the timt-time and space- 
space components of the tensor field W . The evolution 
equation for a(t) in massive gravity has the form 




8(C + 2) 



b 3 -(l-4C)- + §(y 2 -3)]= o;u , n 



y a 



-2fc(fci+4fc 2 )/(fci + k 2 ), C 



(4) 



where a 2 = Aki, ff 2 

(3 2 /a 2 . The perfect fluid form of is assumed and the 
equation of state is taken in the simple form, p — top with 
oj constant. In contrast to GR the function b(t) is not 
arbitrary and gets determined through a(t) itself ; there 
is an algebraic relation between a(t) and b(t) pH] . This 
relation acquires a simple form for £ = 1/4 



(-1 + V7 + 9a 4 ) 
3a 2 



(5) 



II. PHASE SPACE ANALYSIS 

In this section we investigate the phase space behavior 
of cosmological evolution [p9| in finite range gravity for 
both the gravitons behaving like tachyons, i.e, a 2 , 1 < 0. 
For simplicity, we will be working with £ = 1/4. The 
evolution equation ( 4 ) in this case for a >> 1 and 
uj = 0, where interesting things are expected to happen, 
simplifies to 

a' 2 (T) = ^j- + A ™« 2 + ^A m a~ 2 + y^A m a~ 4 - A m (6) 



and for a << 1 



12 A «i -4 i a 

a a + — 



(7) 



where A r 



( HM i 

VS(C+2) J 8tG Po 

derivative with respect to r = y/ (8irGpo)t 
obtain the expression ( 6 ) we have used the following 
series expansion for y(a » 1) 



and prime denotes the 
In order to 



higher order terms 



Before proceeding further, a comment about Eq.( 6 ) is 
in order. From the point of view of GR , Equation ( 6 ) 
can be thought as an effective Friedmann equation where 
the first term on the right hand side of this equation is 
the usual term contributed by the matter density and 
the rest of the terms containing A m originate from the 
Mass tensor appearing in the finite range gravity the- 
ory. The second term behaves like cosmological constant 
and is responsible for drawing the late time accelerated 
expansion in Babak-Grishchuk theory. It is remarkable 
that such a term automatically appears at late stages of 
evolution in massive gravity. The terms containing A m 
smoothly vanish in the mass-zero limit leaving behind the 
usual Freedmann equation. It is surprising that an arbi- 
trarily small mass of gravitons can change the character 
of evolution, turning deceleration into acceleration. The 
third term in Eq. 6 mimics radiation whereas the last 
term looks like curvature term. It is interesting to notice 
that Eq. 6 looks formally similar to the effective Fried- 
mann equation in brane- world cosmology p0| ]. It should 
however be kept in mind that unlike the FRW cosmology 
where the numerical value of the scale factor does not 
carry any physical significance and the scale factor can 
be made larger or smaller than one at any given epoch, 
such a choice no longer exists in massive gravity. As 
already mentioned, the function b(t) is not arbitrary in 
Babak-Grishchuk theory, it gets determined through a(t) 
via Eq. 4 which clearly shows that the numerical value of 
scale factor is indeed important. Eq. 6 is valid for large 
values of " a" and in this regime out of all the terms con- 
taining A m , the second term is most dominant which is to 
compete with the first term allowing ultimately the tran- 
sition to acceleration. It is further important to mention, 
as pointed out by Babak and Grishchuk|l9|, that for no 
values of parameters, a = in the present case ( how- 
ever, for a 2 ,/? 2 > 0, the scale factor goes through regular 
maximum and minimum). At early times, the evolution 
of universe is described by Eq. 7. The first term on the 
right hand side (RHS) of this equation dominates at early 
epochs and mimics stiff equation of state. 

It would be suggestive to recast the expression ( 6 ) to 
look like Newton equation in one dimension fell 



a (r) 



da 



V(a) 



with the first integral of motion 



E 



+ V{a) = 



(8) 



(9) 
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FIG. 1: Plot of the potential V(a) given by Eq.10 for 
A m = 10" 12 . This is a generic form of potential which in 
case of specially flat universe allows deceleration to change 
into acceleration. 



where V(a) is given by the expression, 



V(a) 



A 



6 



rn 2 

a 



7 

12" 



A„,a 



7 . _ 4 A m 

36 Ama + T" 
(10) 

The potential V(a) is plotted in Fig.l. The system, 
in the present case, has sufficient kinetic energy to sur- 
mount the potential barrier and, consequently, the scale 
factor passes through a point of inflection at which the 
deceleration changes into acceleration. In order to ana- 
lyze the dynamics of the system we recast the equation ( 
8 ) in the following form 

v = —uv (11) 



*' = A m - lA m v 4 - lA m v e - *L - u 2 (12) 
6 9 6 



where 



1 a' 

v = - , u = — 
a a 



The constraint equation ( 9 ) takes the form 



u A m A m 2.7. 4,7. 6 v . . 

y = — -— v+ u AmV + m AmV + y (13) 

To find out the critical points we set right-hand side 
(RHS) of Eqs. (11) and (12) to zero and we obtain two 
fixed points 



Ucr = 0, 



vo 



and 



v cr = 0, u cr = +A 1//2 
where Vg is a solution of the algebraic equation 

.,6 i 7 A „,4 



7 7 v 

-A m vl + -A m v* + y - A„ 







(14) 



FIG. 2: Phase portrait of cosmological evolution in massive 
gravity for a 2 ,/3 2 < , i.e., for both the gravitons in the 
theory behaving like tachyons. Trajectories, corresponding to 
different values of parameter A m , starting anywhere in the 
phase space, consistent with the constraint given by Eq. (13), 
end up at the stable critical point (0,1). 



v cr = —(Am) 1 / 2 is not relevant here. At the first critical 
point, a'(r) = and this is not possible in the present 
case i .e, for a 2 , (3 2 < 0. Indeed, omitting terms of higher 
order than v 3 in equation ( 14 ), one gets the approximate 
value of — (SAm) 1 / 3 . This critical point violates the 
constraint equation ( 13 ) and therefore should be dis- 
carded. The second critical point which satisfies the con- 
straint equation ( 13 ) is very interesting, u cr = +(A m ) 1 / 2 
means that 



a'{T)/a(T) = A]l 2 or a(r) oc e 



and v cr = means that the exponential behavior is ap- 
proached asymptotically for large values of a(r). Thus, 
this scenario describes a rapidly accelerating universe 
eternally at late times and like any other generic model of 
quintessence, the finite-range theory of gravity will also 
be faced with the problem of future event horizon fl32[. 
In order to check the stability of the critical point we 
perturb about fixed point (0, +A 1 / 2 ) 



v = v c , 

U = U r 



Sv 
- Su 



Plugging these in 
terms we get 



11 ) and ( 12 ) and keeping only linear 



Sv = 



Su' = —2u rr 5u 



(15) 



(16) 



with the solutions 



Sv 



(x e 



Su 



,-2A x / 2 7 



which clearly demonstrate the stability of the critical 
point under small perturbations. As the system evolves, 



4 



the point moves in the (u, v) plane along the curve given 
by the constraint equation ( 13 ) on which lies the criti- 
cal point. Stability guarantees that the motion will take 
place towards the critical pont. We have evolved the sys- 
tem numerically for different values of the parameter A m , 
we find that trajectories starting anywhere in the phase 
space, consistent with the constraint, end up at the stable 
fixed point as shown in Fig. 2. Therefore the accelerated 
expansion is a late time attractor of finite range gravity. 

To summarize, we have considered the cosmological 
evolution in the massive gravity theory proposed by 
Babak and Grishchuk. We have investigated the phase 
space behavior of the system in case both the massive 
gravitons in the theory are tachyons. We find that late 
time acceleration is an attractor of the system We there- 
fore conclude that the late time accelerated expansion is 



a generic solution of Babak-Grishchuk theory which does 
not require any speculative form of matter to achieve 
it. It would be interesting to confront the predictions of 
Babak-Grishchuk theory with observations and investi- 
gate the constraints imposed by supernova results on the 
massive gravity. 
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